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Abstract 

Motivated by the excess events that have recently been found near the endpoints of the 
double beta decay spectra of several elements, we re-examine models in which double beta 
decay can proceed through the neutrinoless emission of massless Nambu-Goldstone bosons 
(majorons). Noting that models proposed to date for this process must fine-tune either 
a scalar mass or a VEV to be less than 10 keV, we introduce a new kind of majoron 
which avoids this difficulty by carrying lepton number L = —2. We analyze in detail the 
requirements that models of both the conventional and our new type must satisfy if they are 
to account for the observed excess events. We find: (1) the electron sum-energy spectrum 
can be used to distinguish the two classes of models from one another; (2) the decay rate for 
the new models depends on different nuclear matrix elements than for ordinary majorons; 
and (3) all models require a (pseudo) Dirac neutrino, having a mass of a several hundred 
MeV, which mixes with v P . 
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1. Introduction and Summary 



Recently, a mysterious excess of high-energy electrons has been seen in the electron 
spectrum for the double-beta (/3/3) decay of several elements. This kind of observation was 
first made in 1987 for the decay 76 Ge — > 76 Se +2e~ by Avignone et al. [1], although the 
effect was discounted when they, as well as other groups, subsequently excluded a signal 
having the original strength [2]. The mysterious events are back, however, with the UC 
Irvine group now finding excess numbers of electrons near but below the endpoints for 
100 Mo, 82 Se and 150 Nd, with a statistical significance of 5cr [3]. Such events also persist in 
the Ge data [4], [5], at approximately a tenth of the original rate. 

Since these are difficult experiments, it is possible that the anomalous events will turn 
out to be due to systematic error, or to a hitherto unsuspected nuclear physics effect. 
But they may also be the fingerprint of a new fundamental interaction [6], [7], in which 
two nucleons decay with the emission of two electrons and a light scalar, the majoron, 1 
rather than the usual two-electron, two- neutrino decay [8], [9]. If so, these observations are 
of vital importance since they provide us with a glimpse of physics beyond the standard 
electroweak theory. 

We assume for the sake of argument that the excess high-energy electrons are really 
due to majoron emission, denoted by (3(3 m . Our goal is to explore the implications of (3/3 m 
taken together with the other known constraints on neutrino physics. In so doing, we have 
found that the candidate models capable of describing majoron emission from nuclei fall 
into two broad classes. 

In the first class of models for (3(3 m — which to our knowledge includes everything 
that has been proposed until recently [7], [10] - the majoron is the Nambu-Goldstone 
boson associated with the spontaneous breaking of a U(l) lepton number symmetry. The 
only way to get an observable rate in this context is to have either a scalar mass or vacuum 
expectation value (VEV) of the order of 10 keV. We refer to these as 'ordinary' majoron 
models (OMM's), and denote their associated beta decay by f3f3 om . We provide here a first 
comprehensive analysis of which OMM's can give a large enough rate of f3f3 orn . 

In addition, we have recently proposed [11] a second, qualitatively different, sort 
of majoron that does not require such a small scale. Unlike OM's, this new majoron 
carries a classically unbroken lepton number charge, and is the Nambu-Goldstone boson 

1 The term "majoron" was originally used for the Nambu-Goldstone boson associated with spontaneous 
breaking of lepton number, since the same lepton number breaking induced a Majorana mass for the 
neutrinos. We enlarge the meaning of the name in this paper by applying it even if the scalar is massive 
or if the model in question does not generate Majorana masses. 
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for a symmetry distinct from lepton number. We accordingly call such theories 'charged' 
majoron models (CMM's), and denote the associated decay by f3f3 cm . 2 

Our main results, briefly summarized in ref. [11], are: 

• 1: The two classes of models predict different electron spectra for majoron-emitting 
double beta decay, which may therefore be used to identify the type of process that is 
being observed. 3 

• 2: If the majorons are electroweak singlets and the couplings are renormalizable, then 
(3f3 m is observable only if there is a neutrino which mixes appreciably with the electron 
neutrino and whose mass is at least ~ (50 — 100) MeV. CMM's are further constrained to 
have the mass of this neutrino also not much heavier than a few hundred MeV. 

We start, in the following section, with a brief summary of the experimental situation, 
parametrizing the size of the effect in terms of the strength of a hypothetical Yukawa 
coupling between the majoron and the electron neutrino. There follows a formulation of 
the naturalness problem faced by OMM's. This motivates the introduction and definition 
of our alternative: the charged majoron. 

Section (3) proceeds with an analysis of the /3/3-decay rate for a theory with generic 
neutrino masses and Yukawa couplings. We derive the shape of the predicted electron sum- 
energy spectrum for all of the models of interest, as well as a general momentum-space 
parameterization of the relevant nuclear matrix elements as a sum of six form factors. 
General formulae for the (3(3 decay rates are presented in terms of these form factors, which 
we also translate into nuclear matrix elements in the nonrelativistic impulse approximation 
for the weak interaction currents. 

Sections (4) and (5) then apply the general expressions derived in Section (3) to 
specific models of the ordinary and charged majoron type. The properties of the particle 
spectrum required for a sufficiently large (3(3 rate are determined, and the necessity of a 
neutrino with mass M > 100 MeV is explained. We show that, for CMM's, M must also 
not be much heavier than this scale if the observed anomalous (3(3 rate is to be accounted 
for. A similar conclusion follows on less robust grounds from naturalness considerations 
for OMM's. 

Having established what conditions are necessary for producing the observed decay 
rate, we turn in Section (6) to a discussion of the constraints these theories must satisfy to 

A variation on this theme in which this broken symmetry is gauged has been discussed in ref. [12]. 

o 

A similar spectrum can arise for OMM's if two majorons are emitted simultaneously, as in the models 
of ref. [13]. 
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avoid conflict with other experiments. Searches for heavy neutrinos in the decays K — > ev 
and 7r — > ez/ currently furnish the most restrictive laboratory limits. Nucleosynthesis is 
given particular attention in this section, since it would rule out the existence of light 
scalars that are required in both ordinary and charged theories. We show how these 
bounds can be evaded by somewhat complicating the various models. 

2. General Considerations 

We begin by parametrizing the size of the anomalous effect in the data, and expound- 
ing the theoretical naturalness issue which provides the biggest challenge in accounting for 
the excess events. For the purposes of model building, the salient features of the anoma- 
lous events are that they are above where standard ftfav decays contribute appreciably, 
yet below the endpoint for the decays. These facts preclude their interpretation as either 
(5(52v or the neutrinoless j3f3Q V . 

Another crucial input comes from e + e~ annihilation at LEP. The precise measurement 
of the Z-boson width for decay into invisible particles constrains its couplings to putative 
light scalars. Any model in which the rate for Z — > (light scalars) is appreciable, for 
example that of Gelmini and Roncadelli [7], is ruled out. We therefore focus on scalars 
that are electroweak singlets [6] . Although it is possible for majorons to be an admixture of 
both singlets and fields carrying electroweak charges, they have no advantages over purely 
singlet majorons with respect to the beta decay anomalies. In fact these models suffer 
even more severely from the naturalness problems outlined later in this section, and so we 
will not consider them further. 

2.1) The Size of the Effect 

There are currently four experiments measuring double beta decay with sufficient 
precision to potentially see the excess events observed by the UC Irvine group. Two report 
no excess, with one of these quoting an upper bound [14] that is marginally in conflict 
with the Irvine result. 

To compare the effect in various nuclei, we follow the experimental practice of quan- 
tifying the (3(3 m rate using a hypothetical direct Yukawa coupling between the electron 
neutrino and a massless scalar, ip. The rate for majoron emitting double beta decay fol- 
lows from the Feynman graph of Fig. (1), evaluated using the effective interaction 

£phen = ~9eS Ib^e <P- (!) 
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Table 1 lists the coupling strength, g e ^, needed to produce the observed signals in the 
various double beta decay experiments. Our analysis used the nuclear matrix elements 
(M = ([cr n • a m — (g v /g A ) 2 ]h(r))) found in Staudt et al. [15] to estimate the rates for the 
two-neutrino and majoron decay modes. The details of how these matrix elements arise are 
explained more fully in later sections. Here g v and g A are the axial and vector couplings of 
the weak currents to the nucleon and h{r) is a neutrino potential. To quantify the number 
of excess events, we choose (by eye) a threshold energy, .Eth, above which the anomalous 
events begin and the contribution from ordinary (5(3 decay is negligible. The data are taken 
from ref. [3] for the elements 82 Se, 100 Mo and 150 Nd, and from the published spectrum of 
76 Ge in ref. [4]. In all of these cases the excess events comprise R = 2 to 3 % of the total 
number observed. Interestingly, g e ^ lies in the range 8 x 10 -5 to 4 x 10 -4 for all elements. 

Although the coupling apparently needed for 100 Mo looks disturbingly large compared 
to the others, this may be due to uncertainties in the evaluation of the nuclear matrix 
elements. As described in ref. [15], the Ou matrix element for 100 Mo in particular suffers 
from the near collapse of the random phase approximation that was employed. 

We also quote here, for comparison, the results of the Heidelberg-Moscow-Gran Sasso 
group, who claim a 90% c.l. upper bound for 76 Ge of g e g < 1.8-10 -4 in [14]. A similar bound 
of g e ff < 2.0 ■ 10~ 4 is reported for Xe decay by the Neuchatel-SIN-Caltech collaboration 
[16]. 



Element ^(y" 1 ) R E th 9eS 



76 Ge 2 x 10" 23 0.02 1.5 1 x 10" 4 

82 Se 2 x 10" 22 0.03 2.2 8 x 10" 5 

100 Mo 3 x 10" 21 0.03 1.9 4 x 10" 4 

150 Nd 3 x 10" 20 0.02 2.2 2 x 10" 4 



Table 1: The parameters required for emission of ordinary majorons in double beta decay. is the 

inverse half-life of the anomalous events; and R is the ratio of anomalous to the total number of events. 
E th (MeV) denotes our choice for the threshold value of the sum of the electron energies, above which 
essentially only excess events appear. g e ff is the phenomenological coupling (defined in eq. ((1)) required 
to explain the excess rate. 

Besides these laboratory experiments in which the electron energy spectrum is di- 
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rectly measured, there are also several geo/radio chemical experiments. In these the final 
abundance of daughter products is measured, so only the total decay rate can be deter- 
mined. Since the energy spectrum is unknown, it is impossible to directly determine which 
process is responsible for the decay. For comparison with Table 1, we show in Table 2 the 
couplings, <7 e ff, of eq. (1) that would be allowed assuming the total decay rate were due 
to the majoron-emitting process. Since these values are comparable with those in Table 
1, confirmation of the laboratory excess events would likely imply a significant role for 
majoron emission in the geophysical observations. 

The predictions for Te are of particular interest because of a recent measurement 
of the ratio of decay rates C = r( 130 Te) /V ( 128 Te) = (2.41 ± 0.06) x 10 3 . 4 Taking the 
ratio of the lifetimes is useful because some of the uncertainties in their experimental 
determination are expected to cancel. As we will explain in more detail in subsequent 
sections, the significance of this ratio lies in its strong dependence on the relative phase 
space for the two decays [18]. It is therefore sensitive to the integrated electron spectrum, 
which can discriminate between the different possible decay processes. 

The allowed coupling for 238 U is included here for completeness, although we have 
been informed that the discrepancy between the 238 U observations of Turkevich et al. 
[19] and the calculations of Staudt et al. [15] have now been resolved by improving the 
theoretical estimates. 



Element 







9eS 


128 Te 


1 x 10" 25 


0.23 


3 x 10" 5 


130 Te 


4 x 10" 22 


30 


6 x 10" 5 


238-q 


5 x 10" 22 


33 


2 x 10" 4 



Table 2: The parameters consistent with emission of ordinary majorons in double beta decay. T .\ is the 

1 1 2. 

total inverse half-life, assumed to consist completely of anomalous events. Q is the total phase space 
available for each decay measured in units of Ul e . As in Table 1, g e ff is the required coupling, defined in 
eq. (1). The changes in Te relative to ref. [11] reflect the new measurements of ref. [17]. 



4 This result of Bernatowicz et al. [17] was used to constrain f3/3 om emission by W. Haxton at Neutrino 
92, Granada, Spain [18]. 
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2.2) The Naturalness Issue 



One of the first puzzles that must be addressed by any theory of the anomalous events 
is how PP m could be seen without evidence for the neutrinoless decay, (3(3q v . If the emit- 
ted scalar is the Nambu-Goldstone boson for spontaneous breaking of lepton number, as 
in OMM's, then /3/3 0u must exist at some level due to the generation of Majorana neu- 
trino masses. We argue that OMM's answer this question by requiring some dimensionful 
parameter in the scalar potential to be of order 10 keV. 

The small scale arises because the same VEV, u, that breaks lepton number in these 
models typically also generates a Majorana mass for the electron neutrino whose size is 

m UeUe ~ g e s u. (2) 

The Majorona mass gives rise (3(3q v decays which would have been seen if it exceeded the 
experimental limit 

m VeUe < 1 eV. (3) 

Together with the inferred coupling strength, g e ff ~ 10~ 4 , this bound implies an upper 
limit for the lepton-number breaking VEV of 

u < 10 keV. (4) 



One might try to avoid such an artificially small scale simply by having no breaking 
at all, u = 0. In this case (3(3^ is completely forbidden by lepton number conservation. 
The question then becomes why the emitted scalar in (3(3 m should be so light. Since the 
experiments resolve events within 100 keV of the endpoint, the scalar must be no heavier 
than 100 keV. Spontaneous breaking of lepton number naturally satisfies this constraint 
since the Majoron is an exactly massless Nambu-Goldstone boson, but if lepton number 
is unbroken, the smallness of the mass would seem to require fine tuning of parameters in 
the Lagrangian. 

In either case — by kinematics if lepton number is unbroken, or from eq. (2) if it is 
broken — we are led to a mass scale in the scalar sector of the order of 10 — 100 keV. 
Introducing it by hand is at best repugnant. Naturalness demands that the smallness of this 
new scale, relative to the higgs VEV, for instance, must be stable under renormalization. 
Otherwise we have a new hierarchy problem, which is particularly severe if the light scalars 
carry electroweak quantum numbers, as in the triplet majoron model [7]. In that case, loops 
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involving the electroweak gauge bosons generate contributions to the scalar potential that 
are of order y/ot/Air M w > 100 MeV. 

It has been claimed in ref. [10] that a majoron coupling g e g ~ 10~ 4 is small enough 
to generally allow such a hierarchy below the weak scale to be stable. But in the OMM 
that these authors consider, the effective coupling measured in (3(3 m decay is g e s ~ g9 2 , 
where the mixing angle 6 is bounded by neutrino oscillation and decay experiments to be 
very small. This means that the coupling dominating radiative corrections, g rather than 
<7eff, is not small: g ~ O(l). Therefore the corrections to the small 10 keV scale will tend 
to be at least four orders of magnitude bigger than the scale itself and fine-tuning must be 
invoked. 

We show in Section (4) that the scalar hierarchies in these models can be made sta- 
ble under renormalization by taking advantage of the small couplings and masses within 
the neutrino sector, but only in some corners of parameter space having potentially trou- 
blesome phenomenology. For example, the OMM model of Section (4) points to heavy 
neutrinos in the mass range of several hundred MeV that mix appreciably with v e . Even 
though such models are technically natural, they suffer from the aesthetic problem of 
requiring mysteriously small dimensionless scalar self-couplings, £ < 10~ 14 . 

2.3) Introducing Charged Majoron Models 

The above comparison suggests a third option in which the light scalar mass and the 
absence of the neutrinoless decay, Pfiov, can both be naturally understood. To do so, we 
still assume that the emitted scalar is a Nambu-Goldstone boson in order to insure its small 
mass. The absence of f3j3o u is also guaranteed if the spontaneously broken global symmetry 
is not lepton number, which we assume remains conserved. Thus is completely 

forbidden because it is a AL = 2 process. The majoron-emitting decay is still permitted, 
however, provided that the massless Nambu-Goldstone boson itself carries lepton charge 
L = —2. We dub such particles "charged majorons," and show as one of our main results 
that they lead to qualitatively different features for double beta decay, thus allowing them 
to be distinguished from ordinary majorons. 

3. General Properties of the Double-Beta Decay Rate 

Next we derive expressions for the rates of the various possible double-beta decay 
processes. Although a number of excellent reviews exist [20], [21], detailed formulae are 
presented here for several reasons. Our first goal is to highlight the differences in predictions 
between OMM's and CMM's, since the CMM's have not been considered in earlier work. 
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Secondly we want to isolate the dependence of our results on the nuclear matrix elements, 
since these are the most uncertain factors. For generality, we introduce a form-factor 
parametrization of the decay rate which relies simply on the symmetries of the problem. 
Expressions for these form factors in the familiar nonrelativistic impulse approximation 
are subsequently derived. 

There are essentially two properties of double-beta decay that can be measured or 
computed: the shape of the spectrum as a function of the energies of the two emitted 
electrons, and the overall normalization of this spectrum, which determines the total decay 
rate. Only the second of these quantities depends on the size of the nuclear matrix elements. 

Consider the differential decay rate for the four processes to which the experiments 
are potentially sensitive: /?/?2i/> Pflov, f3f3 orn and f3f3 cm . The amplitudes for the first two 
depend on the Feynman graphs of Fig. (2) or Fig. (3), respectively. Those for the majoron 
emitting processes require instead the evaluation of Fig. (1) using the appropriate majoron 
couplings (more about which later). 

It is convenient to write the resulting rates as 

dr(PP) = {GF °™ 3 ec)4 \A((3(3)\ 2 dn((3(3), (5) 

where G F is the Fermi constant, 9 C the Cabibbo angle, A(f3f3) a nuclear matrix element, 
and dQ(/3/3) the differential phase space for the particular process. The observables are 
taken to be the energies of the two outgoing electrons, (k = 1,2). Deriving explicit 
formulae for A(f3f3) and dfl(f3f3) is the goal of the remainder of this section. 

Eq. (5) shows that the decay rate depends on the nuclear matrix elements only as 
an overall multiplicative constant. The only approximation that must be made to derive 
this form from the graphs of Figs. (1-3) is the neglect of the dependence of these matrix 
elements on the final-state lepton energies and momenta. This is a good approximation for 
the neutrinoless modes in which we are interested because the final leptons (plus majoron) 
can carry at most the endpoint energy, Q ~ (1 — 3) MeV, while the nuclear matrix elements 
are characterized by the nucleon Fermi momentum, p F ~ 100 MeV. Corrections to this 
approximation thus introduce a relative error of order Q/p F ~ a few per cent. 

3.1) The Electron Energy Spectrum 

Consider first the electron energy spectrum, dQ(/3/3) in eq. (5). This factor is deter- 
mined solely by the leptonic part of the appropriate Feynman graph. From Figs. (1-3) it 
is straightforward to find the following results. 
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The PPov decay is essentially two-body since the nucleus is too heavy to carry away 
any appreciable kinetic energy. The electron phase space is 



dft(Ppou) =-^2 6 (Q-^~ £ 2) II Pk£kF(e k ) de k . (6) 

A 1 

Here p k = \p k \ is the magnitude of the electron three- momentum, and Q is the endpoint 
energy for the electron spectrum, determined by the initial and final nuclear energy levels, 
M and M', to be Q = M — M' — 2m e . F(e) is the Fermi function, normalized to unity in 
the limit of vanishing nuclear charge. 

In contrast, the phase space for the other three processes can be written in a similar 
form, 

i 2 
d£l{Wi) = ^ (Q ~ ei - e 2 ) ni II P^kF(e k ) de k . (7) 
n k=i 

Only the spectral index rii differs between fifiiv, ft Pom and (3(3 C m decays, 

n 2u = 5; n cm = 3; n om = 1. (8) 

For f3f3 2 K and l3f3 om these values of n; simply reflect the phase space for the corresponding 
process. But for (3(3 C m there are two extra powers of (Q — e\ —€2) due to the proportionality 
of the leptonic matrix element to the majoron energy, a distinctive and generic feature of 
CMM's that we elucidate in Section (3.5) below. We have assumed that the boson emitted 
in [3[3 orn or /3/3 cm was massless; if it has mass m one must use ((Q — ei — e 2 ) 2 — m 2 ) 1 / 2 in 
place of (Q — ei —62) in eq. (7). 

The difference between ri2v = 5 and n om = 1 has long been recognized as a way for 
experimenters to recognize a possible admixture of these types of decays; they lead to 
differently shaped curves for the differential rate, dT/de, as a function of the sum of the 
electron energies, e = e± + €2- The surprising fact that charged majorons have an index 
n cM = 3 intermediate between ftfov and (3(3 om therefore makes it possible, in principle, to 
determine whether a distortion in the ftfov spectrum is due to ordinary or charged majoron 
emission. Fig. (4) shows the shape of the sum-energy spectra for the three possible values 
of the spectral index. 

The spectral shape can also have implications for the total decay rate which, be- 
ing an integral over the sum energy spectrum, depends strongly on n^. Roughly speak- 
ing, each successive power of (Q — e) in dT/de suppresses the total rate by an additional 
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power of (5/(100 MeV). Therefore geophysically-determined decay rates, such as the ratio 
C = r( 130 Te)/r( 128 Te) defined in the previous section, may ultimately prove useful for 
distinguishing between different models. Once the relative strength of Pfev to (3(3 m decays 
is better determined, a definite prediction for Q will become possible. If, for example, the 
decay rate is dominated by the majoron emitting process, then (3(3 om decay predicts too 
small a ratio [18]; we find that (((3f3 m) = 93. This number includes a factor of 5/7 due 
to the ratio of nuclear matrix elements as computed by ref. [15], and the more significant 
factor of (30.4/0.23), due to the difference in phase space for the two decays (see Table 2). 
Because of the small endpoint energy for 128 Te compared to that of 130 Te, the same ratio 
for (3(3 cm decay is much larger: C(PPcm) = 770, and is closer to the experimental value. 

3.2) The Nuclear Form Factors 

The other observable constraining models of majoron-emitting double-beta decays 
is the total rate for any given decay. This requires a knowledge of the matrix element 
denoted A([3[3) in eq. (5), forcing us to deal with the uncertainties in calculating nuclear 
transition amplitudes. The latter can be written as a sum of six form factors, with which 
we parametrize the dependence on nuclear physics. The form factors can subsequently 
be expressed (as we do below) within the context of a given nuclear model. We start 
by defining the form factors, and then use them to specify A{j3j3) for the various decay 
processes in Sections (3.3) through (3.6) below. 

The nuclear matrix element that appears in the evaluation of Figs. (1-3) is 

W a ^P\p)^(2nf ] l^ J J d*x(N'\T*{Ux)M0)]\N)e^. (9) 

Here = ^7^(1 + 75)^ is the weak charged current that causes transitions from neutrons 
to protons, and | N) and | N') represent the initial and final + nuclei in the decay. E and M 
are the energy and mass of the initial nucleus, N, while E' and M' are the corresponding 
properties for the final nucleus, N' . The prefactor, a/ EE' /MM', is required to ensure 
that W a p transform as a tensor since, as is common in the literature, we use nuclear states 
which are not covariantly normalized: (p|p ; ) = 5 3 (p — p'). The (27r) 3 is conventional, and 
is required in order to put our matrix elements into the standard form once the overall 
centre-of-mass motion of the nucleus is separated out. 

A priori the tensor W a p is a function of the four-momenta, and P'^ of the initial 
and final nuclei, as well as four-momentum transfer between the two currents, p^. This 
dependence can be significantly simplified, however. For (3(3qu and (3(3 m , p^ is of the order 
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of the nuclear Fermi momentum p F ~ 100 MeV, whereas the difference (P — P')^ is only 
a few MeV and may therefore be neglected compared to p M . Then the dependence of W a p 
on P M and P'^ may be replaced with the single variable w M , the common four- velocity of 
the initial and final nuclei. For (3 feu, the momentum transfer, p M , is itself also of order the 
energy released in the decay, and so in this case W a p may be simplified even further by 
approximating p^ pa 0. 

It is also straightforward to show that the Bose statistics of the weak currents, J a , 
imply that W a p{u,p) = Wp a (u, —p). Using the aforementioned approximation, the most 
general possible form for W a p is [ll]: 5 

W af3 (u, p) = Wt 7] a p + W 2 U a U{3 + W 3 p a pf3 + W 4 (p a Up + ppU a ) 

+ W 5 {p a Up - Pf3U a ) + iw 6 e a t3apU a p p , (10) 

where the six Lorentz-invariant form factors, w a = w a (u ■ p,p 2 ), are functions of the two 
independent invariants that can be constructed from p^ and it M . Under the reflection 
p — > —p, all the Wi are even except for w 4 , which is odd. 

By evaluating the leptonic parts of the (3(3 matrix elements and contracting with W a p, 
one can show that, to leading order in lepton energies, only its trace, W a a , enters into the 
rates for ftftiv, Pflov, and (3(3 om . In terms of the Gamow- Teller and Fermi nuclear form 
factors — which we define in the nuclear rest frame by w F = Woo an d w GT = £V Wa - 
we therefore retrieve the familiar linear combination 

W a a = w F -w GT . (11) 

For (3f32v we may to a good approximation neglect p M . This permits two important simpli- 
fications: (i) we may drop all but the form factors w\ and W2, and (ii) we may approximate 
these two form factors by constants, Wi(u ■ p,p 2 ) ~ Wi(0, 0). In this limit there is a direct 
relation between w\ and W2 with w F and w GT , given by w\ ~ \w GT and W2 ~ w F +^w GT . 
For (3f3ov and (3(3 om however, p M is large and so W3 and W4 may also contribute significantly 
to w GT and w F . 

The next step is to express A((3(3), and hence the double-beta decay rate dT (/?/?), in 
terms of the form factors, w a . Before doing so, we pause to present explicit expressions for 
these form factors, modeling the nuclear decay as the independent decay of its constituent 
nonrelativistic nucleons. Besides giving some intuition as to the potential sizes to be 
expected for these form factors, these expressions allow a connection between our form- 
factor analysis and the nuclear matrix elements that appear in the literature. 

5 The reader should be advised that we define our form factors here differently than in ref. [11]. 
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3.3) The Form Factors in the Nonrelativistic Impulse Approximation 



The common practice in the literature is to provide expressions for the double-beta 
decay rates with the nuclear matrix elements computed using explicit models of the nu- 
cleus. In this section we present expressions for the form factors using such a model. This 
gives a point of contact between the formalism we present here and the rest of the litera- 
ture. Besides providing a check on our calculations, the expressions we obtain give some 
indication of the size that might be expected for each of the form factors. 

Our evaluation starts by inserting a complete set of states, into the matrix 

element of eq. (9). Working in the nuclear rest frame ( where EE' /MM' = 1 ) and writing 
out the time ordering, we have: 



Contact with the literature can be made once we perform the following approximations: 

• (1) The Closure Approximation: In this approximation a sum over intermediate states of 
the form J2 X F(E X )\X)(X\ is simplified by replacing the X-dependent prefactor, F{E X ), 
by F(E) where E is the energy averaged over the states that contribute to the matrix 
element in question. In the present example we may also use the information that M — M' 
is much less than M and E to replace M' with M throughout. 

• (2) The Nonrelativistic Impulse Approximation: The next simplification is to model the 
nuclear decay in terms of the independent decay of its constituent nucleons, which are 
taken to be nonrelativistic. We work in the position representation, as is conventional 
in nuclear physics. In this representation, the weak currents acting on the constituent 
nucleons takes the following form: 



W aP = (2tt) 3 / d 4 xJ2e ipx (N'\J a (x)\X){X\M0)\N) 0(x ) 



x 



+ (N'\Jp(0)\X)(X\J a (x)\N) e(-x )] 




{N'\Jt3(0)\X)(X\J a (x)\N} 
po - E x + M - ie 



(12) 



J (x)=^5(x 



rn)r+ (g v - g A C n ) + Q(v 2 /c 2 ) 



n 



J(x)=^5(x 



r n)r+ {g A a n - g v T) n ) + 0(v 2 /c 2 ), 



(13) 



n 
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where we have included terms up to 0(v/c) in the nucleon velocities. The sum here runs 
over the constituent nucleons, with the position of the 'nth' nucleon denoted by r n . The 
operator a n similarly denotes the Pauli spin matrices acting on the nth nucleon spin, while 
r+ is the isospin raising operator for this nucleon. As in Section (2.1), g v ~ 1 and g A ~ 1.25 
represent the usual vector and axial couplings of the nucleon to the weak currents. 

The operators C n and D n represent the 0(v/c) contributions to the weak currents, 
and are included here since some of the form factors vanish in the limit that v = 0. They 
are defined in terms of the initial and final four-momenta of the decaying nucleon, (E n , P n ) 
and (E' n , P' n ), the Pauli spin-matrices, cr n , the mass of the pion, m^, and the mass of the 
proton, M p , by [20], 

C n = (P„ + P' n ) • a n /(2M p ) - (E n - K)(P n - Pn) • Vn/ml, 

D n = [(p n + p;) + i^(p n - p' n ) x ay /(2M P ). (14) 

Here up = \{g p — g n ) — 4.7 is a combination of the proton and neutron spin g-factors that 
originates from the contribution of 'weak magnetism.' 

The final step is to separate the overall motion of the nucleon centre-of-mass, R, out 
of the nuclear wavefunction. For a nucleus labelled by its overall momentum, P, as well 
as its other quantum numbers, a, we write: 

giPR 

(ri,...,r A |P,a) = ^ 3/2 (h, ■ ■ ■ , r A \a)), (15) 

where the 'reduced' coordinates, r n , are subject to the constraint J2 n *n = 0. 
These approximations give the following results for w F and w GT : 



2i(ig 

w 



2 



1 U ' nm 



w 



ct - -5 f2 ,. ((N'\ J2 e ~ iP ' rnm ^n r m a n -a m \N}}, (16) 



where \x = E 1 — M is the average excitation energy of the intermediate nuclear state, and 
r nm is the separation in position between the two decaying nucleons. We neglect the 0(v/c) 
corrections to this expression. 

The only other combination of form factors which arise for /?/? 2 ^, flflov, f3f3 om and 
P 'P cm decays are w<s and wq, and these arise only in f3j3 cm . In the impulse approximation 
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we are using, these expressions vanish at lowest order in v/c, forcing us to go to the next 
higher order. We find that 



w * = | | 2 , 2 Z ^V ■ s • «iV|e- ip - r - [g 2 AC n a m - C m a n ) + ^(D m - D n )] |AT»; 

^9a9vP 
\p\ 2 (pl - \x 2 + ie) 



w e = , rJ^^T-, ■ ^ ■ ((N'\e-*- r - [D n x a m + a n x D m ] |iV», (17) 



in which X^mn T m T n are implicit. 

We may now complete the calculation by expressing the various double-beta decay 
amplitudes, A((3fii), in terms of the nuclear form factors w\ through wq. For this purpose 
we must specify the form for the interactions and neutrino masses to be used in evaluating 
Figs. (1-3). We consider each of the four decay processes separately in the following 
sections. 

3.4) The I3l3 2 u Rate 

For completeness we start with the standard two-neutrino decay, fifliv- Evaluating 
the total rate using the leptonic part shown in Fig. (2), and comparing with eq. (5), we 
deduce that the nuclear part of the amplitude is approximately 

~— ^= [4^(0,0) -w 2 (0,0)], 

« -4= MO, 0) - w GT (0, 0)] . (18) 

7TV15 

For simplicity all final lepton energies and masses have been ignored. Thus only the 
form factors evaluated at zero argument appear because, for Pfav decay, conservation of 
momentum and energy determines the nucleon recoil four- momentum, p M , in terms of the 
energy-momentum of the final-state leptons. 

3.5) The (3(3o„ Rate 

To evaluate Fig. (3) for the (3(3q u decay rate, one must know the neutrino mass spec- 
trum. We consider a general mass matrix for an arbitrary set of Majorana neutrinos, 

£mass = -\ Vi {mijiL + mijin) v jy (19) 
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where rriij = rriji is the left-handed neutrino mass matrix, and ^ L (7^) are the usual 
projectors onto left-handed (right-handed) spinors. The physical masses rrii are given by 
the square roots of the eigenvalues of the matrix m^m — not necessarily by the eigenvalues 
of m itself, which may be complex. The electron-flavor row of the associated 'Kobayashi- 
Maskawa-type' matrix for the weak charged-current interactions is denoted by V^. 

With this choice the /3/3ou decay matrix element becomes 

A WM = £ VI mi J **, (y^-) ■ (20) 

Although the range of integration runs over all possible neutrino four- momenta, p M , the 
nuclear form factors w a act to cut the integrals off at the Fermi momentum and energy, 
p F and E F . The contributions from heavy neutrinos thus become suppressed, decoupling 
as 1/m, asm starts to exceed this scale. 

Using the approximations of Section (3.3) for the form factors in W a a leads to the 
familiar Gamow- Teller and Fermi expressions, 

(^4 ( p2 _ m 2 + z£ j = h(r nm ; m t ) (g> - g\a n • a m ) |iV», (21) 

where h(r nm ;mi) is the neutrino potential function defined by: 

r n m;m)=^Jdp u{u + fi) ; u = (p + m ) . (22) 

Again, J2mn T m T t i s implicit in these expressions. 

An important special case is that in which the neutrino masses are negligible compared 
to the nuclear scale, p F . Then one can use the massless propagator in eq. (20) and make the 
replacement rrii = m UeUe , since the integral is to a good approximation independent 

of i. Thus the rate vanishes in the absence of a direct Majorana mass for the electron 
neutrino, as it should. 

3.6) The Rate for Ordinary Majoron Emission 

For majoron-emitting decays we wish to evaluate Fig. (1), and this requires a knowl- 
edge of the neutrino-majoron coupling. For generality's sake we take the form 

£ V w = -\ + bitfn) Vj if* + c.c. (23) 
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If the scalar field is real then (23) still applies, but with the restriction that fry = a*j. 
For example, the phenomeno logical interaction of eq. (1) represents the case of a single 
neutrino with b* VeVe = a VeVe = -ig e s- 

Evaluating Fig. (1) using this interaction and neglecting, as before, the final state 
lepton energies and momenta leads to the amplitude 



I / 



d A p 



W a a {a i:j mimj + p 2 bij) 



(p 2 — mf + ie)(p 2 — m 2 - + is) 



(24) 



For neutrino masses that are much smaller than p F this expression simplifies to the form 



A{PPom) = -4V^ 



VeiVejbij 



f d 4 p f 
J W V 



W a 
p 2 + is 



(25) 



which involves the same combination of nuclear matrix elements as appears in eq. (20) 
for A(f3f3Q V ) with light neutrinos, a result first pointed out in ref. [9]. The sum over mass 
eigenstates simply gives the coupling in the flavor basis, b UeVe , which must vanish in a 
renormalizable theory if the majoron comes from an electroweak singlet field. Thus in 
renormalizable singlet- majoron models it is necessary for at least one neutrino to have a 
mass m > p F ~ (50 — 100) MeV. Such a neutrino can generate an effective, nonrenormal- 
izable b VeV& coupling upon being integrated out. Notice that this observation rules out the 
simplest singlet-majoron model [6], (in which the standard model is supplemented by a sin- 
glet scalar field without additional intermediate-mass neutrino species) as an explanation 
for the anomalous (3(3 events. 



3. 7) The Rate for Charged Majoron Emission 



For charged majorons the interaction (23) must be further constrained to reflect the 
fact that (p now carries lepton number. Suppose that the global symmetry for which cp 
is the Nambu-Goldstone boson acts on the neutrino fields in the following way: 8v = 
i>(<llL — Q T 1r) v i with generator represented by the matrix q. Then, as shown in Appendix 
A, the majoron coupling matrix to neutrinos may be written as 

a = -j(q T m + mq), 

b = +Uqm* +m*q T ), (26) 
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where / is the decay constant, proportional to the symmetry-breaking scale. Note that 
Nambu-Goldstone bosons carrying an unbroken charge are associated with nonhermitian 
generators (for example, the longitudinal component of the W ± bosons), so that q T ^ q* 
in what follows. 

Eqs. (26) are equivalent to the statement that it is possible to redefine the neutrino 
fields in such a way as to ensure that the neutrino-boson coupling has the derivative form 

£<puv = y ^7 M (<m - Q t 1r)v d^ip + c.c. (27) 

The equivalence of this interaction with the Yukawa formulation is demonstrated explicitly 
for double-beta decay in Appendix B. 

The big surprise now comes when eqs. (26) are substituted into the result (24) for 
the (3(3 orn decay rate. As is shown by brute force using the Yukawa couplings in Appendix 
C, the result vanishes identically! This is a reflection of the general statement that the 
amplitude, A{j3j3 cm ), vanishes as the energy of the emitted majoron goes to zero [recall 
that we ignored all final state momenta in deriving (24)], a fact which is most easily seen 
using the variables for which the neutrino-majoron coupling takes its derivative form as in 
eq. (27). 

This result depends crucially on having the emitted Nambu-Goldstone boson carry an 
unbroken quantum number — in this case lepton number. The same result does not apply 
to ordinary majorons, even if they are true Nambu-Goldstone bosons rather than being 
massive. This statement may be puzzling on reflection, since in this case also one can put 
the majoron-neutrino coupling into the derivative form of eq. (27). The resolution of the 
paradox is that for OMM's the rest of the amplitude is singular in the limit of vanishing 
majoron energy, leaving a nonzero result. For the details of this argument we refer the 
reader to Appendix D. 

The upshot is that in CMM's one must work to next higher order in the final lepton 
energies than was done to get eq. (24). The extra factors of the majoron momentum can 
be put into dQ((3f3 crn ), and account for the difference between n cm and n om in eq. (8). In 
the rest frame of the decaying nucleus, the nuclear matrix element turns out to be 6 



ij J 



d 4 p 



p 2 (w 5 + Wq) 



(p 2 — m 2 + ie)(p 2 — m 2 + is) 



(28) 



6 We thank C. Carone for pointing out an error (corrected here) in this equation as it appeared in ref. 
[11]. 
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Whereas previously it was the trace of W a p that arose in the decay rate, here it is the 
skew-symmetric part, parameterized by the form factors w$ and wq, that appear. In 
the nuclear rest frame these form factors are given by = Pi(Woi — Wio)/(2|p| 2 ) and 
w G = e ij kPiW j k/(2\p\ 2 ). 

Using the approximations of Section (3.3) for these form factors leads to the formulae 
of ref. [II]. 7 We note that the neutrino potential that results from doing the momentum 
integral in eq. (28) does not give the usual expression, eq. (22), because of the different 
momentum dependence of the form factors (17). 

In fact, there are a number of important differences between the charged majoron 
amplitude (28) and the corresponding result for ordinary majorons, eq. (24): 

• 1: Eq. (28) depends on completely different form factors than the corresponding expres- 
sion for any other kind of majoron-emitting double beta decay. In fact, we know of no 
variety of (3(3 which depends on 105, 8 and this matrix element therefore appears not to 
have been computed by anyone yet. On the other hand wq would appear in (3(3o u if there 
were right-handed currents, and it has been calculated in ref. [22]. An interesting feature 
of this computation is that the value of wq does not appear to be suppressed by the nucleon 
velocity, v/c, as would have naively been expected, but is instead rather large. 

• 2: In the limit where all neutrinos are much lighter than p F , the flavor dependence of 
the amplitude becomes proportional to the same combination of couplings as appeared for 
OMM's: Yliij VeiVejbij = b UeVe . This direct coupling to the electron neutrino must vanish 
in any renormalizable CMM's, because we have assumed that the Nambu-Goldstone bosons 
to all be electroweak singlets. Analogously to OMM's, it follows that in any CMM at least 
one of the neutrinos must have an appreciable mass: > p F ~ (50 — 100) MeV. 

• 3: As may be seen from eq. (28), if either of the neutrinos in the graph are large compared 
to p F , then the result becomes suppressed by at least two powers of the heavy neutrino 
mass: 1/m 2 . Notice that this is a stronger suppression than the 1/m behavior that follows 
from eq. (24), for OMM's. 



These two matrix elements correspond to what was called A 1 in ref. [11]. We have corrected the 
erroneous coefficient of 7/9 which multiplies A\ there. There is a p-wave contribution to the amplitude 
which we called A 2 , omitted here because it is expected to be much smaller. The amplitude A 3 is also 
omitted here because it can be seen to vanish identically. 

Q 

except for a highly subdominant Coulomb/recoil correction to /3/3 „ that would give an S— P 1 / 2 final 
state for the electrons-see eq. (C.2.12b) of reference [20] 
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4. Ordinary Majoron Models 



We now construct a viable alternative to the original singlet and triplet majoron 
models, since these are not able to yield an observable rate of double beta decay while 
still satisfying all other experimental constraints. In this section we focus on ordinary 
majorons, reminding the reader that here 'majoron' means any light scalar with couplings 
to neutrinos, regardless of whether it is a Nambu-Goldstone boson. It will be shown that 
if PPom occurs at the rate suggested by present experiments, one can infer that the masses 
and mixing angles of the neutrinos are in a range where they are potentially observable by 
other kinds of experiments. 



4.1) A Minimal Model 



It is easy to invent a minimal ordinary majoron theory encompassing the low-energy ef- 
fects of more complicated physics at the electroweak scale. Let ip be a complex electroweak- 
singlet scalar carrying —2 units of lepton number. The lowest dimensional operator cou- 
pling <p to leptons, while respecting the gauge symmetries and global lepton symmetry, 
is 

C = -^(LH)(H T L c )p + c.c. (29) 

Here L = i^f) L and H are the usual left-handed-lepton and Higgs doublets. This inter- 
action can be derived from a more fundamental theory, such as the one given in the next 
section, by integrating out heavy particles of mass M. k is a dimensionless number that 
depends on the coupling constants of the underlying theory. 

Once the Higgs doublet is replaced by its expectation value, (H) = v = 174 GeV, 
eq. (29) reduces to a coupling of the form of eq. (23) , with strength 



a v.v. — 0, b v v — 



KV 



(30) 



The imaginary part of ip therefore couples axially as in eq. (1), with g e g = kv 2 /\/2M 2 . 
Because of the requirement g e ff ~ 10~ 4 , it follows that M/y/n ~ 10 TeV, consistent with 
the assumption that the particles of mass M can be integrated out when analysing double 
beta decay. 

If the light scalar, <p, should also develop a VEV, then the effective coupling of eq. (29) 
also induces an majorana electron-neutrino mass, m v& = g e s (<p), which is consistent with 
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the present upper bound only if (if) < 10 keV. The simplest assumption is that ((f) = 0. 
This illustrates the general arguments of section (2.2) in the present example. 

The electron spectrum for (3(3 om that would be predicted by this effective coupling 
can be consistent with the excess events that are seen, provided that at least one scalar 
mass eigenstate is lighter than 100 keV. This conclusion holds regardless of whether ((f) is 
strictly zero or not, since the decay rate found using the scalar couplings of eq. (29) in the 
general expression of eq. (24) is sufficiently large even for massless neutrinos, as would be 
implied by a vanishing VEV. 

The alert reader may wonder how the above model can lead to observable scalar 
emission even when the neutrinos are massless, since this is in apparent contradiction 
to the general result for (3(3 om decay that was stated in Section (3.5) above. There we 
claimed that f3/3 orn is suppressed if all neutrinos are much lighter than the scale, p F , of 
the nuclear matrix elements. The contradiction is only apparent, however, because the 
argument of Section (3.5) presupposed only dimension-four (i.e. renormalizable) Yukawa 
couplings, and so does not include those of eq. (29). In fact, this effective coupling can be 
obtained by integrating out the heavy neutrino that is required by the general arguments 
in a renormalizable theory, as we demonstrate shortly. 

An imperative question in this scenario is why the potential for <f should contain 
such a small scalar mass or vacuum expectation value. But somewhat surprisingly, the 
hierarchy between this small scale and the weak scale is technically natural in the sense of 
being stable against renormalization, at least within the low-energy effective theory below 
the heavy scale, M. Quantitatively, there are two types of dangerous terms within the 
scalar potential of the effective theory, 

pVV, Z<P*<pH^H, (31) 

whose coefficients must be extremely small, p < 10 keV and £ < 10~ 14 , if m v is to be kept 
< 10 keV. If we choose to define the running of these couplings within the decoupling- 
subtraction renormalization scheme, 9 then both couplings run logarithmically, except for 
the discontinuous quadratic contributions when a particle is integrated out at its threshold. 
The initial conditions for the renormalization-group (RG) equations in this scheme are 
given by the values of the couplings, e.g. p(M) and £(M), at the heavy-physics scale 
\i = M, where the effective theory is matched onto the underlying theory. Provided that 
these initial values are small, the logarithmic RG evolution through the scales p < M in the 

9 This scheme consists of the usual MS scheme, supplemented by the explicit integrating out of any 
heavy particles as the renormalization point is reduced below the corresponding thresholds [23] . 
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effective theory keeps them small. The same is true for the nonlogarithmic contributions 
that arise when the W and Z bosons are integrated out, since these particles couple only 
very weakly to cp. Furthermore, even though the coupling, g e g, to light neutrinos is not 
particularly small, the effects of the operator (29) in loop diagrams are suppressed within 
the effective theory by the small (or vanishing) v e mass. 

Although the small scalar mass is stable within the effective theory below the scale 
H = M, the difficult issue is whether there exists a model for the physics at ji = M, 
which can produce the effective coupling, eq. (29), and still not generate large scalar self 
couplings. Such a question can only be addressed within the context of the underlying 
renormalizable interactions, which are the subject of the next section. 

4.2) A Renormalizable Model 

It is useful to look for a "fundamental" theory whose low-energy limit is the phe- 
nomenological model in the previous section. One would like to know whether such a 
theory exists, whether it has any additional observable consequences, and how much fine- 
tuning it requires. Naturally we seek a candidate with the smallest number of new particles. 
With hindsight, the simplest choice appears to be the addition of a Dirac neutrino, whose 
mass will turn out to be in the range of p F ~ (50 — 100) MeV. Of course we also must 
include the singlet scalar that is emitted in f3(3 decay. The Dirac neutrino can be described 
as two singlet left-handed neutrinos s±, whose lepton number charges are ±1, and the sin- 
glet scalar field must have lepton charge —2. The most general renormalizable couplings 
of the new particles, consistent with the assumed symmetries, are 

C = -\LH^ R s_ - Ms + ^ R s_ - ±g + s + ~f R s + (p - ^g_S-'j R s_ip* + c.c. (32) 

For simplicity we assume that lepton number is not spontaneously broken: (ip) = 0. The 
spectrum then contains three massless neutrinos, v' e , z/ M and z/ r , together with a massive 
Dirac neutrino, v^. The relation between the left-handed weak-interaction eigenstates and 
the left-handed mass eigenstates is 

z/ e = v' e cos 9 + i>h sin 6*, 

*- = v t 

s + = —is' e sin 9 + Uh cos 9, (33) 

where tan6> = Xv/M and v% is the charge conjugate of v h . For v h masses, M h , in the range 
of present interest the universality of leptonic weak interactions requires that 9 < 0.1, very 
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conservatively; thus we have the hierarchy M/Xv > 10, and M h = \/M 2 + X 2 v 2 ~ M. If 
i>h is very heavy compared to p F it may be integrated out, resulting in an effective coupling 
of the form of eq. (29), with k/M 2 = X 2 g + /M 2 . 

There are two light scalars which can be emitted in double-beta decay in this model, 
corresponding to the real and imaginary parts of the complex field, <p. The total rate for 
PPom decay is given by eq. (25), where the Yukawa couplings of the scalar to the light 
neutrino are 

a K< = 0, b KK = sin 2 6g + . (34) 

This latter coupling is also equal to g e g, so the experimentally suggested value of g e s = 
10~ 4 may be obtained by varying the parameters of the renormalizable model in the range 
0.01 < g + < 1, and 0.01 < sin6> ~ Xv/M < 0.1. Setting A = 1, we get an upper limit of 
Mh ~ 10 TeV for the heavy neutrino mass. But if A = 10~ 4 , for example, then Mh ~ 100 
MeV, which is the smallest it can be before the amplitude starts to become suppressed 
by powers of Mh/p F - In that case we must use the more exact expression for the leptonic 
part of the matrix element. This is accomplished by making the replacement 

g+ cos 6 sin 1 



p 2 + ie + \p 2 + ie p 2 - M 2 + is (p 2 - M 2 + ie) : 

M 2 

+ g_ sin 2 #— -h _ (35) 

(p 2 - M 2 + IE) 2 

in eq. (25), which agrees with b v >y e cos 2 9 in the limit of large M h . 
4.3) Naturalness 

Having specified the particle content at the intermediate mass scale, we can now return 
to the question of how natural is the smallness of the scalar masses. We saw in the previous 
section that if the initial values £(M) and p 2 (M) are small, then and p{jj) remain 
small as \i runs to lower energies, for which the effective lagrangian (32) is valid. But this 
by itself is not enough; in addition we must establish whether the matching conditions at 
the heavy-neutrino threshold fx = M are consistent with small values for ^(p) and p(fi) 
at scales p > M. We regard the parameters as being naturally small only if no delicate 
cancellation is needed between their values above (i = M and the quadratic contribution 
arising when Vh is removed from the effective theory. 

The contribution to p(M) and £(M) due to integrating out the heavy neutrino is easily 
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estimated from the graphs of Fig. (5), in which u h is the virtual particle: 



^(M)~^M 2 , 

^ (M) ~iJ (36) 



Contrary to the prejudice that a p = 10 keV scalar mass requires extreme fine tuning, we 
see that it is possible to keep Sp < p = 10 keV and yet have g e ff = 10~ 4 using plausible 
values of the underlying couplings. For example g + ~ 10~ 2 , g_ ~ A ~ 10~ 4 and sin6> ~ 0.1 
implies M h ~ 100 MeV, which then implies, from eqs. (36), 5p 2 (M) ~ (10 keV) 2 and 
6£(M) ~ 10" 18 . 

In the above scenario we have v e mixing strongly with a neutrino in the mass range 
of several hundred MeV; a choice with potentially strong phenomeno logical consequences 
(see Section (6)). However Mh can be pushed to much higher values by letting g_ become 
smaller, since this has no effect on the effective Majoron coupling, eq. (34), that is relevant 
to /3/3 decay. While it might seem less pleasing aesthetically to have g_ <C g + , there are 
no logical grounds for excluding this possibility. In this case the model is safe from any of 
the constraints to be discussed in Section (6) that follow from an MeV scale neutrino. 



5. Charged Majoron Models 



We now repeat the above exercise for charged majoron models, i.e., to construct the 
simplest example both as a low energy effective theory and a renormalizable one. It will be 
seen that our CMM has some close similarities to the OMM just constructed. In contrast 
however we will find that a heavy neutrino in the 100 MeV mass range is not merely 
suggested, but required, in order to achieve a high enough (3(3 crn rate. 

To motivate the specific example, we start with some general considerations. Consider 
the spontaneous breaking of a global symmetry group G down to a subgroup H. The 
resulting Nambu-Goldstone bosons can carry quantum numbers with respect to unbroken 
charges in H only if the original group G is nonabelian, and the unbroken charges do 
not all commute with the broken generators of G. In the Standard Model itself, a global 
nonabelian symmetry acting upon the leptons is precluded by their Yukawa couplings to 
the standard Higgs, or equivalently by the charged-lepton masses. We must extend the 
low-energy particle content in order to devise such a symmetry. In so doing, it is prudent 
to let the new particles be electroweak singlets lest dangerous couplings arise between the 
massless Nambu-Goldstone bosons and charged leptons or electroweak bosons. 
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If the new electroweak singlet neutrinos are integrated out, we obtain a low-energy 
effective coupling of the charged majoron to light neutrinos. It is instructive to write 
down the lowest-dimension such interaction that is possible since this reveals many of the 
features that are common to all underlying models. For CM's the general result that 
the (5(5 amplitude must be proportional to the CM momentum (see Appendix B) suggests 
using field variables for which the derivative couplings are explicit. Because the CM carries 
lepton number L = —2, the usual interaction of the form Velhl^eQ^ is not allowed; the 
current has L = 0. Rather, we need an even number of gamma matrices, 



These are the lowest dimension operators that are possible; note that they are suppressed 
by two more powers of the heavy neutrino mass M than are the OMM effective couplings. 
It follows that A{(5(5 cm ) is suppressed by at least the factor 6 2 qp/M 2 relative to the corre- 
sponding OMM result, where q and p are respectively the average momenta of the majoron 
and virtual neutrino. 

This estimate gives us constraints on the parameters needed if the underlying heavy- 
neutrino model is to reproduce the observed anomalous (5(5 events. Recall the OMM result, 
<7eff(OMM) ~ g0 2 , where g and 9 respectively measure the couplings between the heavy 
neutrino and the majoron, and its mixing with z/ e . Roughly, the corresponding CMM 
result is 



Using Q ~ 1 MeV and p F ~ 100 MeV, and 9 < 0.1, we see that c/ e ff ~ 10 4 , as required, 
only if (i) g ~ 1, and (ii) p F /M ~ 0(1). 

5.1) A Renormalizable Model 

The above considerations may be simply illustrated within a renormalizable model. 
We must first choose the global symmetry group G that will break to give a majoron 
carrying a (7(1) charge. The standard model itself provides us with an example, since if 
SU L (2) x U Y (1) were a global rather than a gauged symmetry, the Nambu-Goldstone bosons 
eaten by W + and W~ would each carry a unit of electric charge. We are therefore led to 
try an analogous global symmetry SU F (2) xU L '(l), which is to be broken down to ordinary 
lepton (electron) number U L (1) by scalar fields (pi. These are like the two components of 
the standard model Higgs in being a doublet under the new SU F (2) symmetry; however 




(37) 




(38) 
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they are gauge singlets. In further analogy, (pi carries a unit of the U L > charge, just as the 
standard-model Higgs carries weak hypercharge. Our field content is completed by adding 
an SU F (2) doublet of right-handed gauge-singlet neutrinos, N±, and two sterile SU F {2)- 
singlet neutrinos s± carrying only the new U L > quantum number, namely L' = ±1. The 
U L /(1) factor is required to permit lepton number to be embedded into the flavour group 
through the mixing of v e and the new singlets. Explicitly, the transformation properties 
of the new fields under SU F {2) x U L >(1) are 

1rN ^ (iV~) ~ (2 '° ); 7*a±~(l,±l); (^ o ") ~(2,-l), (39) 

where the subscripts denote the corresponding charges under the ultimately unbroken 
lepton number, L = — 2T 3 + L'. 

We construct the most general renormalizable lagrangian respecting all the symme- 
tries. The usual standard model particles are taken to be singlets under SU F (2), and their 
U L '(1) quantum numbers are chosen to coincide with their lepton (or electron) number. 
The new mass terms and Yukawa couplings are 

C = -\LH 1r s_ - Ms +1r s_ - g+ (N^ L s + ) $ - g. {N 1l s_) $ + c.c. (40) 

(A direct mass term for iV is forbidden by the SU F (2) symmetry.) Here $ = zr 2 $* is the 
conjugate SU F (2) doublet, with r 2 the second Pauli matrix acting on flavour indices. The 
scalar potential is chosen to ensure that $ gets a VEV, which we assume has been rotated 
to the form 

<•> = (")• <41) 

This breaks SU F {2) xU L >(l) down to U L (1), with the unbroken electron- type lepton number 
symmetry generated by L. 

The mass matrix resulting from eqs. (40) and (41) yields a massless neutrino, v' e , and 
two heavy Dirac neutrinos, ip±, whose masses can be written as 

M ±=\ \m 2 ± \]M 2 -Ag 2 + u 2 (X 2 v 2 + g 2 u 2 ) 
with M 2 = M 2 + AV + (g 2 + g 2 + )u 2 

In terms of left-handed neutrino fields, v' e , ip + and ifj- carry L = +1, while ip' + and 
carry L = — 1. Only the L = +1 fields mix with the electroweak eigenstate, v e , with a 



(42) 
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mixing matrix given by 

^ e = v' e c e + {ip+s a + ip-c a )so, (43) 

with sg = sin 9, eg = cos 9, etc. denoting mixing angles which are given in terms of model 
parameters by: 

tan 9 = 



9-U 

2M X /X 2 v 2 + g 2 _u 2 ( ' 

tan 2a = 



M 2 - X 2 v 2 + (g 2 - g 2 _)u 2 ' 



We are interested in the couplings b^ controlling the f3/3 cm decay. Strictly speaking 
this is a 5 x 5 matrix since there are five left-handed neutrinos, but because lepton number 
conservation only permits the bij coupling among L = +1 neutrinos, the result can be 
simply expressed in terms of v' e and ip±. In the basis (z/, we have 



■ I SgS a SgC a 

bij = — — I sgs a — cesin2ct c#cos2ct I . (45) 



2 



—sgc a eg cos 2a cgsin2a 



This expression has the property that the (3(3 crn amplitude is zero if either of the couplings 
g+ or g_ vanishes. 

In addition to the charged majorons, there is also a neutral one, <ps, corresponding to 
the diagonal generator of SU F (2). Its Yukawa couplings to neutrinos can be read directly 
from the lagrangian once this is expressed in terms of neutrino mass eigenstates. The 
coupling is 

£zw = -—■ va L Vj<fz + ex., (46) 

where c^- is a 3-by-2 matrix whose rows are labelled by the L = +1 states (is' e , if)-), 
and whose columns are labelled by the L = — 1 states V 7 -)- We find 

gsgsp -g_sgCf3 
Cij = -j= I -g-Cgc a s/3 - g+s a C{3 g_cgc a C{3 - g+s a s/3 I . (47) 

^ \-g-CgS a S[3 + g+C a C[3 g_CgS a Ci3 + g+C a S(3 

j3 here denotes the mixing angle amongst the L = — 1 fields, and is given in terms of the 
model parameters by: 

2g + uM , t . 

tan 2/3 = x9 ^ , = 48 

M 2 + X 2 v 2 + (g 2 - g 2 )u 2 v ; 
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We should remark that the mass terms in the ip±, ip' ± field variables have the form 
TjM±ij)±ip' ± + h.c, leading to matrix propagators 



(V>±V±) (V>±V4)A_ i (i M\ 
(V4^±> (V4V4>/ p 2 -Ml\M i)' 

The usual Dirac mass terms M±ip±ip± and propagators i^±^±) = i/{^> — M±) can be 
recovered by making the transformation 

it) i)(tt) w 

where t(j c denotes the usual charge conjugate field. In these more conventional variables, 
only the chirality projections of tfj± have definite lepton numbers. 

5.2) Naturalness 

The naturalness issues are much less severe in this model than they are for OMM's. 
This is because the unbroken lepton number permits the scale, u, of 5 , L r F (2)-breaking to be 
much higher than O(10 keV) without inducing unacceptably large fifioi, decay. The first 
step is to determine how large this scale can be. We saw earlier in this section that the 
conditions for achieving an acceptable (3(3 cm rate require heavy neutrino masses, M± ~ 
p F ~ 100 MeV, with comparatively large heavy-neutrino scalar couplings, g± ~ 1. We also 
found that gu ~ M if the neutrino-scalar coupling is not to be suppressed by additional 
mixing angles, such as a of our explicit example. Taken together, these conditions imply 
an SU F (2) symmetry-breaking scale, u ~ p F ~ 100 MeV. 

Besides being four orders of magnitude larger than the symmetry-breaking scale that 
is permitted for OMM's, the CMM scalar sector is also more natural for another reason. In 
both cases the largest contribution to the scalar potential comes from loops which involve 
the heavy neutrino, of mass M. This contribution is dangerous for OMM's because this 
neutrino is itself much heavier than the lepton symmetry-breaking scale. The same is not 
true of CMM's, however, because for these models both u and M are of the same size. 
As a result, even though these particles couple with nonnegligible strength, g ~ 1, the 
contributions of heavy-neutrino loops are 



^ ' J ' M 



16tt 2 ' 
9 2 X 2 



(51) 
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5p 2 is clearly acceptably small, since all that is required is bp 2 < M 2 . The majoron-Higgs 
coupling, on the other hand, must satisfy St; < 10~ 6 , which is also easily satisfied given 
the phenomenological constraint that V e i ~ Xv/M ~ 0.1, which implies A ~ 10~ 4 . 

6. Other Bounds 

The couplings of majorons to matter are seen most directly in the (3(3 processes which 
have been the main subject of this work, but there are other constraints which must also 
be considered. The most serious of these are laboratory searches for the mixing of v e with 
a heavy neutrino, which is one of the generic predictions of the models we have discussed 
above. In addition, one must take care that majoron emission from stars or supernovae 
does not cause them to burn out prematurely, nor do majorons make so large a contribution 
to the energy density of the universe that they cause too much helium synthesis or cause 
the Hubble expansion to slow too much. In the following sections we discuss the models 
proposed above with regard to these issues. 

6.1) Laboratory Bounds 

It was argued that in order to get an observable rate of (3(3 om or (3(3 cm events, it 
is necessary to have a neutrino v h in the 100 MeV range which mixes with z/ e . Such a 
neutrino could be inferred from a 'spike' it implies for the positron spectra of the decays 
n + ,K + — > e + Vh if it is lighter than the decaying meson. A survey of the Particle Data 
Book [24] shows that pion decay experiments limit the mixing angle to values 9 < 10~ 3 for 
a neutrino with mass M h ~ 100 MeV. On the other hand, the above analysis indicates that 
the minimum angle needed for observable (3(3 m is approximately 0.1 for charged majorons 
and 0.01 for ordinary majorons. 

The bound on the mixing angle from pion decays is easily evaded by taking Mh > 
so that the decay is kinematically forbidden. Note that experimental constraints on mixing 
coming from searches for the decays of Vh do not apply to our models. These constraints 
assume the visible decay channel Vh — > e + e~u e due to weak interactions, but in the present 
situation the weak process is completely subdominant to decays into majorons, Vh — > z/ e </?, 
which would be undetectable. 

One must therefore look to the decays K — > evh for limits on the mixing angle when 
Mh > 140 MeV. The Particle Data Book lists such constraints only up to a mass of 160 
MeV, so one might be misled into thinking that a modest increase in M h above would 
render large mixing angles safe from being ruled out. Actually there exist stringent results 
from KEK [25] that do not appear in ref. [24]. This experiment also restricts 9 < 10~ 3 
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for masses up to M h = 350 MeV. Such a large mass leads to a large suppression of the 
amplitude for (3(3 cm , although not necessarily for j3j3 om . 

An indirect limit on the coupling of v e to heavy neutrinos also comes from tests for 
universality of the weak interactions of leptons in different families. The most restrictive 
test comes from the comparison of electron and muon charged current couplings in pion 
decays [26]. Suppose that v e had mixing angle 9 to a neutrino with mass M > m^, so that 
its weak couplings were suppressed relative to those of by cos 9 (assuming for simplicity 
that does not mix with anything.) The comparison of theory with experiment shows 
that the ratio of electron to muon couplings measured in meson decays is 

^fAl = cos9 = 0.9970 ± 0.0023. (52) 

Taking the one-sigma lower deviation we get a bound on the mixing angle of 

9 < 0.10, (53) 

which is marginally consistent with having observable (3(3 m in our models. Note that a real 
deviation of (52) from unity would be indirect evidence for the sort of neutrino mixing we 
need. 

A further constraint on the majoron coupling to neutrinos comes from searches for 
the decay it — > ev(p [27]. These yield a comparatively weak limit of g e g < 9 x 10~ 3 . 

6.2) Cosmology and Astrophysics 

Because of the weak coupling of the majoron to matter, one might worry that it could 
have deleterious cosmological effects, such as contributing too much energy density if it is 
massive, interfering with the formation of large scale structure due to its decays. Emission 
of majorons from stars or supernovae might also shorten the lifetimes of either. 

In fact the effective coupling g e ff ~ 10~ 4 of a massive majoron to neutrinos is sufficient 
for avoiding the cosmological problems. The majoron lifetime due to the decay <p — > u e u e 
is 

= lGng-^m- 1 , (54) 

which is 10~ 10 s for m v ~ 10 keV, far less than is required by consideration of the density 
of the universe or galaxy formation. 
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In contrast, the majoron coupling to ordinary matter such as found in stars is too weak 
to do any harm. Since the thermal background of neutrinos in a star is negligible, majorons 
are emitted primarily as Bremsstrahlung from electrons. But for the OMM considered 
above, lepton number conservation prevents a single majoron from being emitted; rather 
they must appear in pairs with zero net lepton number. The effective coupling of two 
majorons to electrons is generated by a loop diagram in which a W boson is exchanged, 
Fig. (6). The resulting effective interaction with electrons can be estimated as 

9 2 m e Gp 2 - (rr\ 

16?r 2 M e ^ e > ( 55 ) 

where 9 is the mixing angle between v e and the heavy neutrino, whose mass does not 

— 1/2 

appear because we have assumed it to be much less than G F ' . The amplitude for ip 
emission proves to be some eight orders of magnitude below the observational limit. We 
expect similar results for charged majorons, which must also be emitted in pairs. But in 
addition we need to check the rate of neutral majoron (</? 3 ) emission in the CMM. The 
coupling of (ps to electrons arises at one loop from W and Z boson exchange. To make an 
estimate we have computed only the latter contribution (the two are numerically equal in 
the singlet majoron model [6]). It is shown in Fig. (7). Using the couplings of eq. (47), 
one can eventually find that the effective interaction has the form 

X2 '"' /(0,a,/?)[(l + 2e)ln(l + e)-2] ^ e 75 e, e = (M 2 + - M 2 )/M 2 (56) 



16V2n 2 u 

where /(6>, a, (3) is a function of the three mixing angles of the model — see Section (5.1) 
for their definitions — and which we here conservatively take to be 0(1). Recall that 
Xv ~ OiMi ~ 10 MeV from the requirement of getting observable (3(3 cm . Using the fact 
that Mi ~ u, we get a coefficient of order 10~ 13 . Comparing with the analysis of ref. [9], 
one sees that this is somewhat below the limit from red giant lifetimes of 10~ 6 times the 
electron Yukawa coupling, or 3 x 10~ 12 . 

Because of the higher temperatures in supernovae, weak interactions are in equilibrium 
and there is a thermal population of neutrinos. A coupling of order g e g = 10 -4 between 
neutrinos and majorons is sufficient for bringing the latter into equilibrium as well [28]. 
Therefore, in contrast to the situation for stars, in supernovae majorons are so strongly 
coupled that they are trapped in the core, and do not significantly deplete the normal 
energy flux, in this case due to neutrinos. This will be made more quantitative in the 
next section, where we examine the equilibration of majorons when the universe was at a 
temperature of 1 — 100 MeV: conditions similar to those in a supernova. 
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If majorons are trapped in supernovae they can have an adverse effect on the bounce 
and subsequent explosion [29]. This has only been studied for triplet majorons, using 
restrictive assumptions about the energy dependence of the cross section for vv — > ipip, so 
that no direct conclusions on the models studied here can be drawn. 

6.3) Nucleosynthesis 

A difficulty not so easily surmounted is that the majorons in our models generally 
change the expansion rate of the universe enough to have increased the predicted abun- 
dance of primordial Helium [30]. We will show how this comes about and suggest some 
possibilities for evading the problem. 

Every scalar degree of freedom in equilibrium at MeV temperatures in the early uni- 
verse is equivalent to 4/7 of a neutrino species in its contribution to the energy density, 
and hence the expansion rate. A complex scalar, as in the OMM we have discussed, would 
thus count as 8/7, and the CMM would give 12/7 because it has a total of three Nambu- 
Goldstone bosons. The current limit on the number of additional neutrino species beyond 
those of the standard model is 0.4 [31]. 

In the OMM's, the dominant means for equilibrating massive majorons is the decay 
(p — > vv and its inverse process. The thermally averaged rate is roughly 

T ~ lO" 2 ^™ 2 /T, (57) 

which comes into equilibrium before a temperature of 1 MeV for all scalar masses greater 
than a few eV, assuming g e g = 10~ 4 . Since it becomes increasingly unnatural to have 
scalars lighter than the 10 keV allowed by the (3(3 experimental anomaly, we expect the 
decays to be in equilibrium for massive majorons. 

Charged majorons will suffer fast decays only if they develop a large enough thermal 
mass. Rather than compute this, we focus on the annihilation process vv — > <p<p*. Using 
the interactions of eq. (45), we estimate the thermally averaged annihilation rate to be of 
order 

T ~ -^6> 4 T 5 M" 4 (58) 

in the limit that T <C M. This is some ten orders of magnitude faster than the expansion 
rate at temperatures of an MeV, assuming masses and mixing angles of 9 ~ 0.1 and 
M ~ 100 MeV. In addition, using the neutral majoron interactions of eq. (47), we find 
that the rate for vv — > ipip$ can be suppressed relative to (58) only by a factor of (T/M) 2 . 
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Thus all three kinds of majorons will be in equilibrium at T ~ 1 MeV, in contradiction to 
the nucleosynthesis bound. 

We would like to point out two ways in which the nucleosynthesis may proceed as 
usual, despite the presence of two or three majoron species. One possibility is that the 
tau neutrino mass is close to its experimental upper bound, in the region of 5 to 30 
MeV. If it decays or annihilates into majorons on time scales faster than 1 s, the time 
when neutrinos decouple, there will be one less species of neutrinos, making room for two 
species of majorons, or three with a weak violation of the bound. 

A second possibility is that some neutrino decays into v e plus tp in such a way as to 
heat the electron neutrinos relative to the other species. It was shown that this occurs if 
Vfj, or v T has the desired decay with a lifetime in the range 6 x 10~ 4 s < r < 2 x 10~ 2 s [32]. 
The overpopulation of v e results in prolonged equilibrium between neutrons and protons, 
which compensates for the extra density of the universe in its effect on helium synthesis. 
This idea can be generalized to the annihilations of sterile neutrinos in the mass range of 
a few MeV as well. In fact it is not necessary that the decaying or annihilating particle go 
directly into z/ e 's; as long as it produces particles that are in equilibrium with u e , after the 
decoupling of neutrinos from electrons, it will accomplish the same thing. 

As an existence proof for these mechanisms, we show how the ordinary majoron model 
of section (4.2) can be generalized to include a heavy tau neutrino. Let there be one 
additional sterile neutrino, S3, whose lepton number is the same as that of s_. When we 
include the other two generations, the straightforward extension of the lagrangian (32) 
yields a mass matrix of the form 
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(59) 



in the basis (v e , v T , s_, s + , S3). It is easy to see that the spectrum consists of two 
massless states which are mostly v e and v^, a Dirac neutrino of mass ~ rrii consisting 
mostly of v T and s_, and a Dirac neutrino of mass ~ Mj which is mostly s + and S3. 
Supposing that the intermediate Dirac mass is of order 10 MeV, for example, we see that 
the constraints on mixing angles can be satisfied: 

9 eT = [cm\ — sm 4 ) / (cms — stuq) < 0.01 (60) 
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from searches for peaks in the it — > ev spectrum [26], and 

9 es = [(cm,4 + smi) — 6 eT (cm§ + sm^)/ M ~ 0.01 (61) 

from the requirement that g e g ~ 10 -4 in /3/3 om . Here s/c = M\jMi is the tangent of 
the mixing angle in the sterile neutrino sector, M = (M-j 2 + Ml) 1 / 2 is approximately the 
mass of the heaviest state, and 6 ex denotes the mixing angle between v e and the mass 
eigenstates that are mostly v T or the heavy sterile neutrino. 

It turns out that the tree level couplings that would cause the decay v T — > v e (p 
vanish; nevertheless the annihilation process v T v T — > ipip* goes at a rate comparable to 
(58), however without the mixing angle suppression. The annihilations are therefore very 
efficient in depleting the v T population, as long as the heavy neutrino mass scale M is 
significantly smaller than 100 GeV. Moreover the resulting majorons are still in equilibrium 
with the light neutrinos, so we have the z/ e -heating mechanism in addition to the elimination 
of v T . We note that a tau neutrino in this mass range would not necessarily have manifested 
itself in supernova 1987a through the delayed signal of its decay products. Because it 
interacts so strongly with majorons, its neutrinosphere will be farther out in the core 
where the temperature is lower and the Boltzmann suppression is greater, contrary to the 
usual case where v T is emitted at a higher temperature. Thus the flux of v T : s would be 
greatly reduced relative to the electron neutrinos. 

7. Conclusions 

Motivated by experiments suggestive of majoron emission in double beta decay, we 
have proposed two kinds of models that are able to account for this effect without lepton 
number violation. In the first proposal, the boson is not of the Nambu-Goldstone variety 
but rather has a small mass, which can nevertheless be natural in a technical sense discussed 
above. The second proposal is to let the majoron be truly massless, but carry lepton 
number charge. Coincidentally, both of these schemes suggest the existence of heavy 
isosinglet neutrinos in the mass range of several hundred MeV with significant mixing to 
the electron neutrino. These heavy neutrinos could manifest themselves in the decays 
K — > ve or by nonuniversality in the weak interactions of electrons versus other leptons. 
The models can be consistent with nucleosynthesis constraints if the tau neutrino is in 
the 1 — 10 MeV mass range, or there exist additional sterile neutrinos with a mass of a 
few MeV. The anomaly in the double beta decay spectra, if confirmed, would thus be the 
precursor to several new phenomena in neutrino physics. 
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Note Added 

Since completing this work we have been informed of evidence that the anomalous 
events reported by the UC Irvine group may be due to resolution problems for the higher- 
energy electrons [33]. 

Appendix A. The Yukawa formulation of NGB Couplings 

In this Appendix we derive the general form for the Yukawa coupling of a Nambu- 
Goldstone boson: eq. (26). 

To this end consider an arbitrary set of Yukawa couplings between a collection of spin 
1/2 and spin particles. Such particles may always be cast as Majorana fermions, and 
real scalar fields, $ a . The most general form for their mutual couplings is 

£yuk = ~ 7 ,V J $a + C.C. (62) 

Suppose also that this lagrangian is invariant with respect to the following global 
symmetry transformations: 



stf = id a ( qa y jlL - (ft,)'/?* 



5<$> a = ie a {Q a ) a % b , (63) 

in which both sets of matrices, q a and Q a , are hermitian, and Q a must also be imaginary. 
Invariance of the Yukawa couplings is expressed by the identity 

(q a ) T T a + T a q a + r b (Q a ) b a = o. (64) 

Any explicit left-handed fermion mass matrix, (mo)ij, must similarly satisfy the relation 
q T m + m q = 0. 
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This symmetry is spontaneously broken when the scalar fields acquire their VEV's, 
v a = ($a)) an d the resulting Nambu-Goldstone boson directions in scalar-field space, (p a , 
are given by the action of the symmetry on v a : 

(5 GS $) a = (F" 1 )^^- (65) 

The real, symmetric normalization matrix, is chosen to ensure that the scalar kinetic 
terms remain properly normalized. That is, <9 M $ a <9 M $ a = <9 M </? a d^(p a H , provided that 

( F -i)«7 [ V T Q^Q XV ] (F- 1 )^ = 6 a(3 . (66) 
The Yukawa coupling for <p a therefore becomes 

C yuk = -^r a 7 ,V (Q a v) a (F-Y P <pp + c.c. 

= {q T a Y a + F a q a ) v a (F _1 ) a/3 ^ + c.a, (67) 

where eq. (64) was used in writing the last line. The expression for the right-handed 
coupling follows simply from taking the complex conjugate of this expression. 

Eq. (67) gives the most general form for Nambu-Goldstone boson couplings. It can 
be recast into the form of eq. (26) using some additional simplifying features of the models 
which we consider. Suppose first that no symmetry-invariant fermion mass terms exist, 
mo = 0. Then the Yukawa coupling matrices, T a , which appear in eq. (67) can be traded 
for a dependence on the fermion mass matrix using 

my = r?- v a . (68) 

Next, suppose that the Nambu-Goldstone bosons carry an unbroken U(l) charge, as is 
the case for CMM's. It is then convenient to work with complex combinations of the <p a : s- 
If, for example cpi and <f2 form a multiplet under the unbroken U(l), then the symmetry 
transformations become diagonal when expressed in terms of cp = (cpi + i(p2)/y/2. The 
same steps as before once more lead to eqs. (67), with the proviso that the corresponding 
broken charge, q = (q± — iq2)/V%, need no longer be hermitian. 

The simplest case is if there is only one Nambu-Goldstone boson with a nonzero 
charge, as in the models we consider. Then the normalization matrix, -F -1 , cannot mix 
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ip with any of the uncharged Nambu-Goldstone bosons, and must be proportional to the 
unit matrix in the charged-scalar sector. Denoting the proportionality constant by: 1/ f = 
(F -1 ) 11 = (F -1 ) 22 , we obtain eqs. (26), as required. 

Appendix B. The Equivalence of Derivative and Yukawa Formulations 

A famous property of Nambu-Goldstone bosons is that they only couple derivatively. 
Here we make this property explicit for the couplings of the Nambu-Goldstone bosons that 
are considered in eq. (26) by showing the equivalence of these two formulations for the 
double-beta decay rate. 

Nambu-Goldstone bosons can only couple derivatively because if these fields are taken 
to be constants they completely drop out of the lagrangian density. This is because the 
Nambu-Goldstone directions in field space are defined by performing a field-dependent 
symmetry tranformations on the vacuum, as in eq. (65). For constant fields these trans- 
formations are really symmetries, and so produce no effect at all in the lagrangian. <p a 
only appears to the extent that it varies in spacetime, and so it must couple only through 
its derivatives. 

To see this in the present case, consider the following field-dependent redefinition of 
the fermion fields: 

# = -i{q alL - ql lR )^ (F" 1 )^ pp. (69) 
The fermion mass term changes by 

5£ m ass = -%^{q T a m + mq a )^ L i) (F _1 ) a/3 <pp + c.c, (70) 

where m = T a v a . Notice that this is exactly what is required to cancel the Yukawa coupling 
of eq. (67). A similar cancellation occurs for all of the nonderivative interactions of (p a . 
It is important to notice in this regard that if the broken symmetry should transform 
other particles like the electron, in addition to neutrinos, then these other particles must 
also participate in the field redefinition, eq. (69), in order to remove all nonderivative 
(/^-dependence. 

The </? a -dependence is not completely eliminated, however, since the fermion ki- 
netic terms are not invariant under a spacetime-dependent transformation such as that 
of eq. (69). It is a simple exercise to show that, under the assumptions leading to eq. (26), 
the variation of the kinetic term is given by eq. (27). The latter has been expressed in a 
way that holds even if the generators, q, are not hermitian, as is appropriate for charged 
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Nambu-Goldstone bosons. Since these two forms for the Nambu-Goldstone boson interac- 
tion are related by a field redefinition, they must give equivalent scattering amplitudes. 

To concretely verify the equivalence of these two expressions for the Nambu-Goldstone 
boson couplings, we compute the double-beta decay rate using both eqs. (26) and (27). 
Although the graph in which the Nambu-Goldstone boson is emitted by the neutrino line, 
as in Fig. (1), is not equivalent by itself, we will show that the result becomes equivalent 
once Fig. (1) is added to the remaining graphs of Fig. (8). We denote the leptonic part 
of the amplitude computed from Figs. (1) and (8) using derivative couplings by M£ u , and 
the same amplitude using Yukawa couplings by M^ v ' . 

First consider the evaluation of these graphs using derivative couplings. For generality 
we work with an arbitrary set of majorona fermions, and real scalar fields, and assume the 
following form for scalar and charged- current interactions which appear in the Feynman 
rules for Figs. (1) and (8): 

C^w = ftl^Vij + A {j 75)^' W» + c.c, (71) 

where, on general grounds, A and A are both symmetric matrices, while V and V are 
antisymmetric. For the charged- current weak interactions we take A = V oc T_, where 
T_ is the SU L (2) lowering operator. For the Nambu-Goldstone boson associated with the 
charge q, we have A = \(q + q T ) and V = \{q — q T ). The corresponding Yukawa couplings 
will be denoted by the matrix A oc cv~f L + by R , with a and b given by eq. (26). 

Up to a common overall normalization the leptonic part of the integrands for the three 
graphs become (in an obvious matrix notation) : 

M£"(l) = u( Pl )Y(V + iA l5 )S( Pl - h)4(V + iA<y 5 )S(k 2 - p 2 )Y(V + iA<y 5 )u c (p 2 ), 
Mg v (8a) = u( Pl )i(V + iA l5 )S( Pl + q)^(V + iA l5 )S(k 2 - p 2 ) 7 "(^* + iA* l5 )u c (p 2 ), 
M^(8b) = u( Pl )Y(V + iA lb )S( Pl - hWiV* + iA* lb )S(- P2 - q)i(V + iA<y 5 )u c ( P2 ), 

(72) 

where S(p) = [ip' + m'y L +m*7 H ] 1 is the fermion propagator, thought of as a matrix in 
Dirac and flavour space, while u (u c ) is the (conjugate) electron spinor. 

These expressions can be related to the Yukawa expressions by applying the following 
easily proven identities: 

S(pi - fci)irf(V + iA l5 )S(k 2 - pa) = S(pi - fci)A(V + iA l5 )S(k 2 - p 2 ) 



38 



+ Sfa - h)(V + M75) - (V + iA lb )S{k 2 - p 2 ), 
u(piMV + iA l5 )Sfa + q) = ufa) [(V - iA l5 ) + AS fa + q)} (73) 
S(-P2 ~ q)H{V + iA~f 5 ) u c (p 2 ) = [-(V + iA-y s ) + S(-P2 - g)A] u c (p 2 ). 



The last two of these identities rely on using the Dirac equation for the initial and final 
spinors, ufa) and u c (p 2 ). 

Using these identities in eqs. (72) relates the derivative-coupling and Yukawa-coupling 
results for each graph: 



M£"(l) = M£"(l) + ufa)^(V + iA l5 )Sfa - h)(V + iA l5 )r(V* + iA* l5 )u c fa) 
- ufa)^{V + zA 75 )(V + iA l5 )S(k 2 - p 2 )Y(V* + iA* lb )u c (p 2 ), 
M£"(8a) = M£"(8a) + ufa)(V - iA lb )^{V + iA l5 )S(k 2 - p 2 )Y(V* + iA* l5 )u c (p 2 ), 
Mg v (8b) = M^(8b) + ufa)<y»(V + iA l5 )Sfa - h)Y(V* + zA* 75 )(-V + iA lb )u c fa). 

(74) 



We see that although the result using the two formulations of the Nambu-Goldstone 
boson couplings do not agree graph by graph, their sum is the same provided that the 
scalar- and charged-current coupling matrices satisfy the following conditions: 



[V, V*] + [A, A*] = [V*, A] + [V, A*] = [V, V] + [A, A] = [V, A] + [A, V] = 0. (75) 

These are trivially satisfied if the charged-current generators commute with the charge 
that is associated with the Nambu-Goldstone boson, as is required by the invariance of 
the charged-current interactions under the spontaneously-broken global symmetry. The 
equivalence of the two formulations for double-beta decay is thus established. 



Appendix C. The CM. Electron Spectrum 



It is argued in the text that the vanishing of the double-beta decay matrix element 
A(PPcm) as the majoron momentum goes to zero is a key feature of charged majoron 
models. Here we demonstrate this property in some detail. 

The vanishing of the amplitude is straightforward when the majoron couplings are 
expressed in derivative form, as in eq. (27). In this case the conservation of electric charge 
and lepton number precludes any derivative coupling between the electron and the charged 
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majoron, in the absence of exotic, electrically-charged fermions. As a result, neither of the 
graphs of Fig. (8) contribute to the CMM double beta decay rate. The remaining graph, 
Fig. (1), manifestly vanishes for zero majoron momentum becaus of the derivative coupling. 

It is more complicated to see this result in the Yukawa-coupling language. Lepton 
number conservation forbids a direct coupling between the charged majoron and the elec- 
tron, so the only graph to be considered is again that of Fig. (1). As might be expected 
from Appendix B, however, the result for this graph need not vanish for zero majoron 
momenta until the contributions from all of the relevant intermediate neutrinos have been 
summed. 

At zero majoron energy, the (3(3 m decay amplitude is given by eq. (24), whose integrand 
is proportional to the v e -v e element of the following matrix in flavor space: 



Integrand oc 



1 



jr — m*m 



(m*am* + p 2 b) 



1 



p z — mm" 



h °° 1 



^^(m*m) k (bmm* — m*am*) (mm*) 1 



k=0 



(76) 



As in previous expressions, m = m T denotes the complex left-handed neutrino mass matrix, 
while a and b are the Yukawa coupling matrices of eq. (23). 

The last expression simplifies drastically once eq. (26) is used, which contains the 
information that the majoron is a Nambu-Goldstone boson. After a pairwise cancellation 
of all but one of the terms in the sum over /c, we find that 



Integrand oc 



qm* + m*q 7 



1 



n=0 Kl/ ' 



2\n+2 



\q(m*m) 



n+l 



m* + (m*m) n+1 m*q T ] 



■ (77) 



The significance of this final result lies in the fact that each term in it is proportional to 
a u e matrix element, q Ue j, of the Nambu-Goldstone boson charge. The final point to be 
established is that, for CMM's, all such matrix elements are zero. We are therefore forced 
to work to next order in the majoron momentum, eq. (28), in order to get a nonvanishing 
contribution. 

In order to see why q Va j must vanish in CMM's, consider the symmetry transforma- 
tions in the basis of weak-interaction eigenstates. Then the invariance of the gauge inter- 
actions under the global symmetry implies that q can only transform the entire doublet, 



40 



(^ e ) , into other doublets having the same hypercharge. But since the Nambu-Goldstone 
boson charge, q, has embedded in it two units of the unbroken lepton number, such a 
transformation cannot be made without introducing exotic isodoublet fermions. 

Appendix D. The O.M. Electron Spectrum 

Here we wish to show that the OMM double beta decay amplitude, unlike that for 
CMM's, is nonvanishing even at zero majoron momentum. This is particularly easy to see 
using the Yukawa form for the Nambu-Goldstone boson couplings, which can be directly 
read off from the lagrangian of a given model. In this form the majoron typically couples 
only to neutrinos, and not to electrons. Thus only the graph of Fig. (1) contributes. In 
contrast to CMM's, it is possible to have q Ve j ^ in OMM's (see the previous Appendix) 
and so the decay rate at zero majoron momentum need not vanish. 

The puzzle is to understand this result when the amplitude is expressed in terms 
of the derivative couplings, since in this formulation all of the graphs of Fig. (1) and 
Figs. (8) are explicitly proportional to the majoron momentum, k. The resolution turns 
out to come from the contributions of Figs. (8). For these graphs, in which the majoron 
is emitted from the electron lines, the internal electron goes on shell in the limit as k — > 0, 
causing a singularity in the propagator. The coefficient of this singularity is proportional 
to the vector part of the electron-majoron coupling. (The same singularity leads to the 
familiar infrared divergence of the analogous photon bremsstrahlung graphs in Quantum 
Electrodynamics.) Consequently the electron propagator behaves as 1/k for small k, which 
cancels the explicit /c-dependence due to the majoron's derivative coupling. 
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Figure Captions 



• Figure (1): The Feynman graph which gives rise to double-beta decay accompanied by 
the emission of a Majoron. The four-fermion vertices are those of the usual charged-current 
weak interactions. 

• Figure (2): The Feynman graph which gives rise to ordinary two-neutrino double-beta 
decay as occurs in the Standard Model. 

• Figure (3): The Feynman graph which gives rise to neutrinoless double-beta decay, with 
no Majoron emission. 

• Figure (4): The number of decay electrons as a function of the sum of the two electrons' 
energy. The solid curve represents the two-neutrino decay, the dotted curve gives the OMM 
decay, and the dashed curve gives the CMM decay. All three curves have been arbitrarily 
assigned the same maximum value for the purposes of comparison. 

• Figure (5): The two most dangerous Feynman graphs contributing to the light scalar 
couplings when the heavy neutrino is integrated out. 

• Figure (6): The Feynman graph through which the effective electron- Majoron interaction 
is induced. 

• Figure (7): The Feynman graph which mixes the Z boson with the 'neutral' Goldstone 
boson. Once the Z is attached to a fermion line this induces an effective electron- Goldstone 
boson interaction. 

• Figure (8): The remaining Feynman graphs which contribute to double-beta decay ac- 
companied by Majoron emission. These graphs only arise if direct electron-Majoron cou- 
plings exist, as is the case for OMM's in the variables for which the Majoron is derivatively 
coupled. 
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